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Abstract. In 1958 Jan Oort remarked that the lack of old clusters in the 
solar neighborhood (SN) implies that clusters are destroyed on a timescale of 
less than a Gyr. This is much shorter than the predicted dissolution time of 
clusters due to stellar evolution and two-body relaxation in the tidal field of the 
Galaxy. So, other (external) effects must play a dominant role in the destruction 
of star clusters in the solar neighborhood. We recalculated the survival time of 
initially bound star clusters in the solar neighborhood taking into account: (1) 
stellar evolution, (2) tidal stripping, (3) perturbations by spiral arms and (4) 
encounters with giant molecular clouds (GMCs). We find that encounters with 
GMCs are the most damaging to clusters. The resulting predicted dissolution 
time of these combined effects, idis = 1.7(Afi/10'* Mq)°'^^ Gyr for clusters in 
the mass range of 10^ < Mi < 10^ Mq, is very similar to the disruption time 
of tdis = 1-3 ± 0.5(^^^/10"* Mo)0-S2 Qyj. ^^Yiat was derived empirically from a 
mass limited sample of clusters in the solar neighborhood within 600 pc. The 
predicted shape of the age distribution of clusters agrees very well with the 
observed one. The comparison between observations and theory implies a surface 
star formation rate (SFR) near the sun of 3.5 10~*° M0yr~*pc~^ for stars in 
bound clusters with an initial mass in the range of 10^ to 3 10^ Mq. This can be 
compared to a total SFR of 7 - 10 xlO"*" M0yr~'^pc~^ derived from embedded 
clusters or 3 — 7 10~^ MQyr~*pc~^ derived from field stars. This implies an 
infant mortality rate of clusters in the solar neighborhood between 50% and 
95%, in agreement with the results of a study of embedded clusters. 



1. Introduction 



Jan Oort pointed out in 1958 that "in general, galactic clusters do not appear 
to grow much older than 0.5 Gyrs" (Oort 1958). Later, Wielen (1971) derived 
a mean dissolution time of 0.2 Gyr from the age distribution of open clusters 
within about 1 kpc from the Sun. Spitzer (1958) suggested that the short lifetime 
may be due to the fact that low density clusters (< IMqpc"^) may be destroyed 
by encounters with passing interstellar clouds. 

The destruction of star clusters occurs in two stages (see Bastian & Gieles, 
these proceedings, for a review.): 

(1) At young ages (< 10 or 20 Myr) most of the clusters dissolve because the 
gas that is left over from the giant molecular cloud (GMC) is removed from the 
cluster by stellar winds and supernovae. This loss of binding energy causes a 
large fraction of the embedded clusters to dissolve into the field. This "infant 
mortality" (Lada &: Lada 2003) depends critically on the star formation efficiency 
and appears to be independent of the initial cluster mass. (Bastian et al. 2005 
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suggest that this occurs within ~ 10 Myr, whereas Fall et al. 2005 suggest that it 
occurs over a much longer timescale of ~ Gyr). Theoretical studies, e.g. Kroupa 
8z Boily (2002) and references therein, show that the gas removal phase and the 
resulting infant mortality last about 10 Myr (see also Lada &: Lada 2003). 
(2) At later ages clusters can be destroyed by external effects such as the tidal 
field, and perturbations due to encounters with GMCs, with spiral arms and the 
Galactic disk. This effect is mass dependent: massive clusters survive longer 
than low mass clusters. 

This paper deals with the second effect, but we also derive the infant mortality 
rate for clusters in the solar neighborhood (SN). 

Starting with the pioneering work of Spitzer on the dynamics of star clus- 
ters, many studies have been devoted to the understanding and prediction of 
the survival times of clusters open and globular clusters in our galaxy. These 
theories predicts that the dissolution time of clusters depends on their initial 
mass: massive clusters survive longer than low mass clusters (e.g. Spitzer 1958; 
Wielen 1985; Chernoff Sz Weinberg 1990; Gnedin & Ostriker 1997 and references 
therein). Baumgardt & Makino (2003) (hereafter BM03) showed from A^-body 
simulations that the predicted dissolution time of a cluster in the tidal field of the 
galaxy depends on the number of stars, N, in the cluster as t^is ~ /3 {N/ \n^N)^, 
where is the Coulomb logarithm and the constants (3 and x depend on the 
initial concentration of the cluster. Gieles et al. (2004) showed that this can be 
written as 



tdis = i4(Mi/lO^M0)T (1) 



with 7 = 0.62, where t4 is the disruption time of a cluster with an initial mass 
Mi = 10^ Mq. This same power-law dependence was derived empirically from a 

study of cluster samples in four galaxies by Boutloukos & Lamers (2003) (here- 
after BL03). The dynamical models of BM03 result in a ^4 = 6.9 Gyr for clusters 
in the SN. Since the mean mass of the open clusters in the SN is about a 10^ 
Mq, the "mean lifetime" of clusters in the SN, predicted by BM03, would be a 
few Gyr, i.e. much longer than the empirical value of 0.2 Gyr derived by Wielen 
(1971). 

To solve this discrepancy, we have started a series of studies of different 
effects that might play a role in limiting cluster lifetimes. These effects are: 

(1) mass loss by stellar evolution, 

(2) dissolution by two-body relaxation in the tidal field of the Galaxy, 

(3) perturbations due to the passage of spiral arms, 

(4) perturbations due to encounters with GMCs. 

At the same time we rederived empirically the cluster dissolution times of open 

clusters in the SN, based on the new cluster sample of Kharchcnko ct al. (2005). 
We compare the predicted and the observed age distributions in clusters to see 
which one of the possible destruction mechanisms is the most important one in 
the SN. This comparison also results in an estimate of the star formation rate 
and infant mortality rate. (For the full study, see Lamers &; Gieles 2006.) 
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Figure 1. Top: The surface density distribution of open clusters in the 
Kharchenko et aL (2005) sample, projected onto the Galactic plane. Bot- 
tom: The ratio between the number of old (> 0.25 Gyr) and young (< 0.25 
Gyr) clusters as function of distance. Both distributions are flat within their 
uncertainties up to 600 pc. 



2. The observed age distribution of clusters in the solar neighbor- 
hood 

Kharchenko et al. (2005) published a catalog of 520 galactic open clusters in 
the SN with the values of angular sizes of cluster cores and coronae, heliocentric 
distances d, E{B — V), mean proper motions, radial velocities and ages. These 
parameters have been determined by homogeneous methods and algorithms in- 
cluding a careful procedure of cluster member selection. The basis of this study 
is the ASCC-2.5 - All-Sky Compiled Catalog of about 2.5 million stars down to 
y ~ 14 (completeness limit at y ~ 11-5), with proper motions and B^V mag- 
nitudes based on the Tycho — 2 data and supplemented with Hipparcos data 
and some ground-based catalogs. Cluster membership is based on a combined 
probability which takes into account kinematic (proper motion), photometric 
and spatial selection criteria. Cluster ages were determined with an isochrone- 
based procedure which provides a uniform age scale. This resulted in the most 
homogeneous catalog of open clusters in the SN. Lamers et al. 2005a (L05) have 
shown that the lower mass limit of the clusters in the Kharchenko et al. sample 
is about 100 Mq. 

Figure [U shows the distance distribution of the density of clusters projected 
onto the Galactic plane, in number per pc^. Within the statistical uncertainty 
the surface density is constant up to at least 600 pc. The lower part of the figure 
shows the ratio between old (t > 0.25 Gyr) and young {t < 0.25 Gyr) clusters 
as a function of distance. Up to a distance of about 1 kpc there is no significant 
change in this ratio within the statistical uncertainty. This is important for our 
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Figure 2. The age histogram in units of number per year, in logarithmic 
age-bins of 0.2 dex, of 114 open clusters within d < 600 pc from Kharchenko 
et al. (2005). In order to show the effect of binning the data, the distributions 
are plotted for two sets of bins, shifted by 0.1 dex, with and without squares 
respectively. The error bars indicate the la statistical uncertainty. The dis- 
tribution decreases to older ages, with a small bump around log(i/yr) ~ 8.6. 
For log(t/yr) < 7.5 the distribution is uncertain due to large error bars. 



study because it shows that the sample of clusters within 600 pc has no age 
bias. Figure [2] shows the age distribution in number per year of the 114 clusters 
within 600 pc in the Kharchenko sample. This distribution decreases with age, 
apart from a small local maximum around log(f/yr) ~ 8.6. The distribution 
at young ages is sensitive to the choice of the age-bins and shows a significant 
scatter. The steep slope at log(t/yr) > 8.8 demonstrates that cluster disruption 
is important. 



3. The evolution of the mass of star clusters 
3.1. Mass loss by stellar evolution 

The stellar mass loss from clusters has been calculated by various groups. We 
adopt the GALEV models for clusters with a Salpeter type mass function in 
the range of 0.15 < M/Mq < 85 (Anders & Fritze-v. Alvensleben 2003). These 
models are based on stellar evolution tracks from the Padova group, which in- 
clude overshooting, mass loss due to stellar winds and supernovae. LOS have 
shown that the fraction of the initial cluster mass that is lost by stellar evolu- 
tion, qev{t) = AM /Mi, can be approximated accurately by 

loggev(i) = (logt - aev)^°^ + Cev for t > 12.5 Myr. (2) 

with t in yrs and Oev = 7.00, bev = 0.255 and Cev = —1.805 for solar metallicity 
models with a Salpeter mass function. 
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This function describes the mass loss of the models at t > 12.5 Myr with an 
accuracy of a few percent. The mass loss at younger ages is negligible because 
stars with > 30 M© hardly contribute to the mass of the cluster due to 
the small number of stars. (These stars do play an important role in the infant 
mortality because they are responsible for the fast removal of the gas from the 
clusters.) The mass loss from clusters by stellar evolution can then be described 
as 



3.2. Mass loss by the galactic tidal field 

BM03 have calculated a grid of A''-body simulations of clusters in circular and 
elliptical orbits in the tidal field of a galaxy for different initial cluster masses, 
galactocentric distances R, and different cluster density profiles. The stars follow 
a Kroupa initial mass function without primordial binaries, and stellar evolution 
is taken into account during the evolution. BM03 and Lamers et al. (2005b) 
have shown that the predicted dissolution time can be expressed as a function 
of the initial cluster mass as 



where ^4 is a constant that depends on the tidal field strength of the galaxy in 
which the cluster moves and on the cllipticity of its orbit. We adopt the value 
of t4 = 6.9 X 10^ yr from BM03 for clusters in circular orbits at Rq = 8.5 kpc. 
Applying this result of BM03 to clusters in the SN, which are on average of 
lower mass than considered in their study, results in an overestimation of the 
dissolution time of low mass clusters by less than a factor 2 or so (see Fig. 5 of 
BM03). This is because low mass clusters, with total lifetimes less than about 
a Gyr, contain massive stars during most of their lifetime. These stars are very 
efficient in ejecting stars by two-body interactions. 

The mass loss due to the Galactic tidal field can then be written as 



3.3. Mass loss by spiral arm perturbations 

When clusters cross a spiral arm, the enhancement of the ambient density gives 
rise to time dependent tidal forces that can accelerate stars out of the cluster. 
These external perturbations are most destructive for clusters with low density 
and for passages with low relative velocity, Vdrift • The solar neighborhood is close 
to the corotation radius of the spiral arms, where the effect is largest. Although 
the time between spiral arm passages at that location is long, oc V^^ift' mass 
loss per passage is high, i.e. oc 

^drift (Gieles et al. 2006a, hereafter GAPZ06). 
GAPZ06 calculated the dissolution time of clusters in the SN due to pertur- 
bations by spiral arms by means of N-body simulations. For the density of the 
clusters, they adopted the observed mean mass-radius relation of clusters as ob- 
served in nearby spiral galaxies by Larsen (2004): rh = 3.75 (M/W^Mq)^ with 




(3) 



tdis = i4 (Mi/lO^Mo)' 



0.62 



(4) 
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A = 0.10 lb 0.03, where is the half mass radius. Here we have converted the 
cluster radius Teff given by Larsen into the halfmass radius using r^s = 0.75rh 
(Spitzcr 1987). For a density contrast of the gas component of the spiral arms 
based on observations, and a mean drift velocity of Vdrift = 12.5 (Dias &i Lepine 
2005), GAPZ06 found a mean dissolution time of 

= ^» (lO^) (^)' = (M/10^Me)'-« Gyr. (6) 

Notice that for A = 0.1 the disruption time due to shocking by spiral arms has 
almost the same mass dependence, i.e. oc M°'^, as that due to the tidal field 
(BM03) and as found empirically (BL03), viz. M^-^^ 

The mass loss of clusters due to spiral arm shocks is 

/dM\ -M(t) / M(t) 

3.4. Mass loss by giant molecular cloud encounters 

Spitzer (1958) suspected already that encounters with GMCs play an important 
role in the destruction of clusters in the SN. Similar to the case of spiral arm 
perturbations, the disruptive effect of the encounters depends on the density of 
the star cluster and on the relative velocity between the cluster and the cloud. 
However in the case of an encounter with a GMC, the mass loss also depends 
on the mass of the GMC and on the distance of the passage at the moment 
of closest approach. This implies of course that the dissolution of clusters due 
to encounters with GMCs can only be described in a statistical way, when all 
effects are properly averaged over the encounter probabilities. 

Gieles et al. (2006b, GPZB06) studied the encounters between GMCs and 
clusters with TV-body simulations. They showed that the relative mass loss due 
to a shock, i.e. AM/M, is only about 20 % of the relative energy gain AE/E. 
This is because a large fraction (~ 80 %) of the energy gained in an encounter 
goes into ejecting stars with a velocity exceeding the escape velocity. This implies 
that the disruption time should not be defined as the timescale for injecting the 
same energy as the cluster binding energy (e.g. Ostriker et al. 1972), but as the 
time scale to bring the cluster mass to zero. 

The disruption time of clusters by encounters with GMCs does not depend 
on the mass of the individual clouds, but only on the mean density of molecular 
gas in the SN. This is because the energy gained by the cluster per individual 
encounter is proportional to the cloud mass Mdoud; but the number of encounters 
per unit time depends on the number density of the clouds. The number density 
times the mass of the clouds is the mean density of the molecular gas. In other 
words: 1000 clouds of 10^ Mq have statistically the same disruptive effect as 10 
clouds of 10^ Mq if we ignore the tidal field (but see §3.5). 

GPZB06 derived an expression for the energy gain and the resulting mass 
loss for the full range of encounter distances, from head-on to distant encounters. 
They adopted a mean midplane density of molecular gas in GMCs near the sun 
of pn = O.O3M0pc"^, a surface density of individual GMCs of En = 17OM0 pc~^ 
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Figure 3. The mass evolution of a cluster with an initial mass of 10^ (left) 
and 10"'' Mq (right) in the solar neighborhood. The mass loss due to the 
four separate effects is indicated. Encounters with GMCs are the dominant 
dissolution effect in the solar neighborhood. A cluster of Mj = 10^ Mq 
dissolves in 0.35 Gyrs and a cluster of 10^ Mq in 8.2 Gyrs. 



(Solomon et al. 1987) and a mean velocity dispersion of clusters and GMCs of 
(j„ ~ 10 kmsec"^. 

With these data they derived a mean dissolution time (tcMc) for clusters by 
GMC encounters in the SN, taking into account the relative velocity distribution, 
the distribution of impact parameters and gravitational focusing. Adopting the 
mean mass radius relation of clusters (Larsen 2004) (see §3.3) they found 



'^^^ = '-^ (lo^k) (^)' = ' Gy^' 



The statistical mass loss rate of a cluster due to encounters with GMCs is 

0.3 

\ dt J GMC *GMC 

Notice that the mass dependence is the same as for dissolution by spiral arm 
shocking, but that the effect is ten times stronger. 



fdM\ -M{t) M{t) _ 1 



3.5. The predicted mass evolution of clusters in the solar neighbor- 
hood 

The mass loss from clusters, due to the combined effects of stellar evolution, 
tidal stripping, spiral arm shocks and encounters with GMCs, is 

dM _ fdM\ 

~dt ~ \dfj + \dfj tidal '^Vdfj ,pi, ^ \dfj GMC ■ 

In writing this expression we assumed that the four effects act independently 

of one another. This may not be the case. For instance, a cluster that has 
encountered a GMC has an increased radius and is thus more susceptible to 
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Figure 4. The predicted dissolution times of clusters in the solar neighbor- 
hood due to the combined effects of stellar evolution, tidal field, spiral arm 
shocks and encounters with GMCs, as a function of the initial mass. Full line: 
total dissolution time. Dashed line: time when the remaining mass is fOO 
Mq. Squares and dotted line: dissolution time due to stellar evolution and 
the Galactic tidal field only, predicted by BM03. Cross with error bar: the 
value of t4 empirically derived by LOS. 



tidal stripping than a cluster that did not have a recent encounter. However, we 
expect this effect to be small, except for the lowest mass clusters. 

We have solved Eq. [10] numerically for clusters of different masses. The 
results are shown in Fig. [3] for clusters with an initial mass of 10^ and 10^ Mq. 
Notice that encounters with GMCs are the dominant dissolution effect in the 
solar neighborhood, contributing as much as the three other effects combined. 
Using the calculated mass decrease of clusters of different initial mass we can 
predict the age of the cluster when it has a remaining mass of 100 Mq, i.e. 
the detection limit for clusters in the SN in the Kharchenko sample (see L05). 
Figure m shows the ages of clusters when their remaining mass is and 100 Mq 
as a function of the initial mass. The almost linear part can be described by 



tdis = 1.7(Mi/10^ Mq)°-6^ Gyr for 3.5 ~ log Mi/ Mq ~ 5 (11) 



The steep decrease at low mass is because clusters with Mj < 10^ Mq quickly 
reach the detection limit of 100 M©. The figure also shows the dissolution times 
by the Galactic tidal field, predicted by BM03 for clusters with an initial con- 
centration factor Wq = 5 in a circular orbit at i?o = 8.5 kpc and the empirical 
value of t4 from LOS. Our predicted timescales are about a factor 5 smaller than 
those due the tidal field only (BM03) but agree with the empirical value (L05). 
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Figure 5. The observed age distribution of an unbiased sample of clusters 
with M > 100 Mq in the solar neighborhood within 600 pc (Karchenko et 
al. 2005; LOS) in units of nr yr~^ is given by squares with the Poisson error 
bars. The full line shows the predicted distribution for a cluster sample with 
a maximum mass of Mmax = 3 10'' Mq and a SFR of 4 10^ Mq Myr"! within 
600 pc from the Sun. The dotted and dashed lines are for Mmax = 1-5 10'* M© 
and Afniax = 6 10"* Mq respectively. The bump around 0.4 Gyr suggests an 
increased SFR at that time. 



4. Comparison with observed age distribution of clusters in the solar 
neighborhood 

Given the initial mass distribution of the clusters, their formation rate, CFR(t), 
and the time it takes for a dissolving cluster to fade below the detection limit, we 
can predict the distribution of observable clusters as function of age or masfl 
Here we are interested in a sample of clusters with M > 10^ Mq. 

For a constant CFR and a power law cluster IMF with a slope of —a = —2 
(Lada & Lada 2003) the number of clusters with M > 100 Mq as function of 
age is 

NM>ioo{t) = C (Mn^(t)-^ - M-l), (12) 

where Mii^{t) is the initial mass of clusters that reach M(t) = 100 A/q at age t. 
Clusters of age t with a smaller initial mass have M < 100 Mq by now. Afmax 
is the maximum initial mass of the clusters that are formed. The constant 
C is related to the star formation rate (SFR) in bound clusters as SFR = 
Cln{Mm^jMu^) for -a = -2. 

Figure [5] shows a comparison between the observed age distribution of clus- 
ters within 600 pc with the predicted distribution for Afmax = 3 10^ Mq. This 
value of Afmax is adopted because the observed distribution shows a steep drop 



^ L05 have derived an expression for the general case of a cluster sample that is set by a 
magnitude limit for any formation history 
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at log i ~ 9.5 (with only one cluster in the last bin) and Fig. shows that this 
corresponds to Mj = 3 10^ Mq. 

The flattening of the predicted distribution at the low age end is due to 
the fact that clusters with an initial mass in the range of about 100 to 300 Mq 
quickly reach 100 Mq (see Fig. U]). The bump in the observed distribution 
around log(t) ~ 8.6 is due to a local starburst (see LOS and Piskunov et al. 
2006). There is good agreement in the shapes of the predicted and observed 
distributions! 

The vertical shift that is applied to the predicted curve to match the ob- 
served one gives a value of C = 10~^'^^ in Eq. [121 which corresponds to a SFR of 
4 10^ Mq Myr"^ for bound clusters in the range of 10^ < M^/Mq < 3 10"^ within 
600 pc from the sun and a surface formation rate of 3.5 lO"^'^ M0yr~^pc~^. We 
can derive the infant mortality rate in the SN by comparing this value with 
the total starformation rate of at least 7—10 lO"^'^ M0yr~^pc~^, derived from 
embedded clusters with a mass 35 Mq, by Lada & Lada (2003) or the rate 
of 3 — 7 10~^ M0yr~^pc~^, derived from the field star population by Miller & 
Scalo (1979). Our value is a factor 2 to 3 smaller than the lower limit of Lada 
&: Lada and a factor 10 to 20 smaller than the one derived by Miller &: Scalo. 
This implies that at least ~70 and possibly even 0.95 % of the stars in the SN 
are born in clusters that dissolve within ~ 10 Myr. 



5. Conclusions and discussion 

Our calculated dissolution times of clusters in the solar neighborhood are about 
a factor five smaller than predicted by BM03 for clusters in the tidal field of our 
Galaxy, with only stellar evolution, binaries and two-body relaxation taken into 
account. This is mainly due to encounters with GMCs. So we can expect that 
clusters in environments with a high density of GMCs will be destroyed very 
effectively. This is especially the case in interacting galaxies, such as M51 and 
the Antennae galaxies. Indeed, Gieles et al. (2005) derived a disruption time 
of only 0.2 Gyr for clusters in M51 at a Galactocentric distance of 1 < i? < 5 
kpc. This is more than a factor 10 smaller than expected on the basis of the 
tidal field only. On the other hand, in galaxies with a small GMC density the 
dissolution time will be set mainly by the tidal field. This is confirmed by the 
comparison between empirical and predicted dissolution times of several galaxies 
by Lamers et al. (2005b) (see Fig. 6 of Bastian &: Gieles, these proceedings). 
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Note added after submission 

In a recent paper Whitmore, Fall & Chandar (astro-ph 0611055) argue that 
our determination of the disruption time of clusters in the solar neighborhood 
is influenced by selection effects. They claim that "not including the young 
embedded clusters from Lada & Lada (2003) can produce an apparent bend in 
the age distribution, which can be misinterpreted in the BL03 models as evidence 
for a specific value of the disruption time" . This is a curious statement because 
the embedded clusters are all younger than 3 Myr and the large majority of 
these clusters will not survive infant mortality (Lada &: Lada 2003). Therefore 
they should not be included in a study of the dissolution of bound dusters. 

Figure O compares our age distribution (open squares with error bars) with 
the one adopted by Whitmore et al. (filled symbols). They adopted an age 
distribution which consists of the sample of embedded clusters of Lada & Lada 
(i.e. with ages < 3 Myr) and open clusters (i.e. which have removed their 
gas) from Battinelh & Capuzzo-Dolceta (1991, MNRAS 249, 76). Lada & Lada 
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Figure 6. Open squares: the age distribution of clusters in the SN (our 
sample). Filled symbols: the age distribution adopted by Whitmore et al. 
(2006). Filled triangle: the embedded clusters of Lada & Lada (2003) at 
log(i/yr) < 6.5. Filled circles: the open clusters by Battinelli & Capuzzo- 
Dolceta (1991). The two distributions agree in the range of overlap. Full 
line: our fit to the data of non-embedded clusters that was used to derive 
the dissolution time of clusters that survived infant mortality. Dashed line: 
Whitmore's linear fit with slope of -1, shifted vertically to match the data at 
log(i/yr) = 7.75. This fit does not agree with the flattening of the observed 
distribution at log(i/yr) < 7.5 and neglects the decrease at log(i/yr) > 9. 



(2003) combined these to show the steep drop (1 dex in log(dA'"/dt) going from 
the embedded to the unembedded phase (i.e. from < 3Myr to > 5Myr), which 
they refer to as "infant mortality" within 10 Myr. 
We see that: 

a. ) the age distribution of the non-embedded clusters from Lada & Lada (2003) 
agrees with the one we derived from the Kharchenko et al. (2005) sample in the 
range 6.75 < log(t/yr) < 8.25. 

b. ) the one high point at log(t/yr) = 6.25 is due to embedded clusters, most of 
which will dissolve by infant mortality (Lada &; Lada 2003). 

c. ) the straight line adopted by Whitmore et al. is a very poor fit to the observed 
age distribution as it seriously overpredicts the distribution at 6.7 < log(t/yr) < 
7.5 and does not explain the drop in the number of clusters at log(t/yr) > 8.0. 
(The peak in the age distribution around log(t/yr) ~ 8.6, probably due to an 
increased star formation episode (L05), is higher than both fits.) 

Therefore, we conclude that the statement by Whitmore et al. (2006), that 
our derived short disruption time is due to the neglect of the embedded clusters, 
is incorrect. In fact, Fig. [6] shows that the linear fit of slope -1, proposed 
by Whitmore et al. in support of their model for a mass independent infant 
mortality over an extended period of time, ~1 Gyr, fits the data of the non- 
embedded clusters poorly. Our fit of the distribution of non-embedded clusters 
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agrees much better with the observations and results in a dissolution time of 
bound clusters that agrees very well with the predictions described in this paper. 



